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Abstract 

The paper contains a proposal for an energy and time represen- 
tation. We construct modes that correspond to fuzzy distributions 
around discrete values of energy or time. The modes form an or- 
thogonal and complete set in the space of square integrable functions. 
Energy and time are self adjoint in the space spanned by the modes. 
The widths of the modes are analyzed as well as their energy-time 
uncertainty relations. The lower uncertainty attainable for the modes 
is shown. We also show times of arrival for massless particles. 



The Pauli theorem revisited 

Two arbitrary states of an elementary system can be transformed into each 
other by symmetry operations. This opens the door to express what can be 
observed of the system, i.e. the system properties, in terms of the generators 
of these symmetries. In the case of the Poincare group they are the momenta 
P M and the angular momenta and boosts M^ v . For the system to be elemen- 
tary they are constrained by the mass shell condition P 2 = m 2 and by the 
spin condition W 2 = m 2 s(s + 1). (The Pauli-Lubanski vector is denned as 

= e fJiUal 3P u M al3 , with % = 1 in this paper unless otherwise specified). 

Notice now the conjunction of both, the four vector character of the 
momenta on one side, with the necessity of introducing conjugate operators 
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to the three-momentum to formulate dynamics on the other. This calls for 
the introduction of a four vector operator conjugate to P^ such that 

ig*" (1) 

Regretfully, this is not possible to attain. The reason is that would 
produce translations in the momentum. If P M is defined on the mass shell, 
then 

P' 1 * = exp (-i5pQ) P» exp (i5pQ) = F" + 5p M (2) 

and, no matter the value of the four vector parameter Sp, P' ^ m . The 
transformations generated by Q pull the momentum P out of the particle 
mass shell. 

The problem is independent of the form of dynamics in use. It has far 
reaching consequences. First, it is necessary to abandon a four vector Q M and 
hence the hope of building a covariant form dynamics in terms of momenta 
and positions p. Then, with time and position demoted to the role of mere 
parameters, it is necessary to introduce quantum fields with creation and an- 
nihilation operators playing the role of conjugate pairs. To our accustomed 
eyes they look like the appropriate recipe for particle number non conserva- 
tion |5j but, Is it so? At least two important lessons remain: 1. That time is 
along the direction that would be conjugate to the solution of the mass shell 
constraint and 2. That time runs due to the constraint. 

With the usual non covariant choice, where the mass shell condition reads 
as P° = H(P), the problem turns into a version of the Pauli theorem [3], 
namely 

H' = exp (—iSET) H exp (iSET) = H + 5E, with5E e R (3) 

This implies that the spectrum of H has to be the real line, something that 
runs again the existence of the physically necessary ground state. So, in 
the time form of dynamics the problem looks like the incompatibility of T 
with the boundedness of the Hamiltonian. On the other hand, if T is not 
selfadjoint, then [H, T] ^ % and there are way outs from the problem. 

Perhaps the simplest case for an observable time is the time of arrival 
(TOA) of a particle at a certain position in one space dimension 0. Classi- 
cally t(x) = m(x — q)/p where q,p are the dynamical variables of phase space 
(initial position and momentum) and x the arrival position. Obtaining the 
quantum version t(x) of this TOA is complicated due to the presence of 1/p 
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and to operator ordering. It is well known that there is no selfadjoint i(x), 
its most symmetrical expression being 



i(x) = — exp(-ixp) J— qJ— exp(ixp) (4) 

V p V p 

The eigenvectors of this operator, \t,x,s) can be given in the momentum 
representation where 1/p is easier to deal with, as 

\p\ P^ 

(P \t, x,s) = 9(sp) \ — exp(z— t) (p \x ) (5) 
V m 2m 

The degeneration parameter s can take the values +1 for right movers, -1 
for left movers. This is the only residue left at one space dimension of the 
continuous manifold of directions present for higher D. 

The bad news come in the form of nonorthogonality of the time eigen- 
states. This is a consequence of the fact that the Hamiltonian is bounded 
from below: a(H) = [E Q , oo) with E = 0. By using (|SJ) we get 

1 f°° 

(t,x, s\t',x, s) = —5 SS > / exp(iE(t — t )) = 

Z7T JE =0 

1 i 1 i 1 

— S ss > lim = - 5 aa > (S(t - V ) (6) 

2tt e->o+t'-t + ie 2 v v ; vr t - V ' y ' 

To get orthogonality it is necessary to move E to — oo something that 
resembles the Weisskopf-Wigner trick for resonances. In any case it is clear 
that the problem source is in a(H) ^ cr(t). 

Marolf devised a procedure, used in ref jH] , to surmount nonorthogonality. 
The idea is to replace t and H by "regularized" expressions that avoid the 
difficulties that arise when p — > 0, namely, to introduce 

Tfl 771 1 p I 

feip) — {"rrif \p\ > e i e ^ se — o~ } withe small and positive (7) 
\P\ e 

Then, 



t(x) = - exp(-ixp) yjfeip) q \f fe(p) exp(ixp) (8) 
an expression tailored for the momentum representation where 

(p \t, x,s) e = 6{sp) -= — exp(iE e (p) t) (p \x) (9) 

'hip) 
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The "regularized energy" is given by 

Ee(p) = £ = {E(P) - E(e) if |p| > e, else E(e) ln(§||)} (10) 

Thus, while -E(p) = when p = 0, E e (p) — > — oo as p — * 0. This solves the 
problem by removing the lower bound in cr{H). It is straightforward to show 
that 

e {t,x,s\t',x,s') e = 5 ss <5(t-t') (11) 

The procedure works fine for the case of free particles. With minor ob- 
vious modifications, it works equally well for relativistic particles and for 
any number of space dimensions. However, it is only suitable for those cases 
where momentum remains constant. So, its very definition brings this pro- 
cedure to a dead end. 

Soon after the publication of ref j3] Giannitrapani jH] observed that the 
time of arrival was an instance of generalized observable endowed with a 
probabilistic interpretation as a positive operator valued (POV) measure. In 
fact, the time of arrival eigenstates form a complete set: 

(p|{£ / dt\t,x,s)(t,x,s\}\p') = (p\p') (12) 

Not being orthogonal they can not constitute a projector valued measure, but 
- as pointed out in ref [U] - nothing prevents from using them to construct a 
POV measure. Giannitrapani showed that the probability that the time of 
arrival at x of the state p be in the range [T, T'] is 

rT' 

dt \t, x, s ) (t, x, s | 



P X ([T,T']) = Tr 



(13) 



The mean value and variance of the TOA are given in |6, , where the problems 
arising from nonorthogonality when trying to obtain kinematical uncertainty 
relations between energy and time are discussed. The paper also analyzes 
the free particle TOA operator of (J1J acting on the domain of infinitely dif- 
ferentiable functions over the compact subsets of values of p G R — {0}. 
In this case the variance of t turns out to be computable as for ordinary 
observables and the Heisenberg uncertainty relation holds. To summarize: 
Endowed with the POV measure interpretation, the TOA became an use- 
ful instrument appropriate to extract physical information from one of the 
primary laboratory events: "when" a detector clicks. No wonder this is the 
customary approach in the current literature. 
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From eigenvectors to wave packets 



Wigner introduced [7] in 1932 a distribution at midway between the position 
and the momentum representation with the aim of describing particle proper- 
ties in phase space. The non positivity of the distribution was an obstruction 
to its use as a probability distribution. Later on, Husimi |8 j introduced sets 
of minimal uncertainty states \qp) in position and momentum: 

(x\qp) = (27HX 2 )- 1 / 4 exp(- ^ ~ ^ ) exptipx) (14) 

The distribution is centered at the point (q, p). Given a system in an arbitrary 
state \ip ), the probability that the system occupy a region in phase space cen- 
tered at (q,p) of half widths (a, l/2c) is given by Ph(q,p) = (2vr) _1 |(gp \ip )| 2 

We now return to the time of arrival of a free particle at a specific position 
in one space dimension. We learnt the difficulties that arise in this seemingly 
simple problem. They can be articulated through the Pauli theorem in two 
complementary ways: The impossibility of finding a conjugate pair of time 
and energy selfadjoint operators or, the difficulty to connect the two different 
spectra a{t) = R and cr(H) = R+. To avoid them, we will follow a procedure 
that resembles Husimi's. 

We start by introducing two functions g u (t) and f v (u) defined over the 
real line t G R = (— oo, +oo) and the positive real line to G R + = (0, +oo) 
respectively. We take v as a positive integer v > 1. Also, we assume g v and 
f v to be Fourier transforms of each other. Finally, 



These functions on which t (respectively uS) act multiplicatively, present 
the following nice property under Fourier transformation: 

. . . dfJu) 

tg v {t) - -i-^ 

. dg v (t) 



dt 



oofv(uj) (16) 



We point out that for v > 1 both derivatives act on these functions as 
selfadjoint operators (/^(O) = 0). It is also remarkable that g u G £ 2 (R) and 
f v G L 2 (R + ). So, they can be given a probabilistic interpretation. Finally, to 
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the canonical pair (t, i ^) acting on t (the would be "Time Representation") 
corresponds the unitarily equivalent pair {—i-^j, acting on u (the would 
be "Energy Representation"). 

This looks like a promising starting point for constructing true conjugate 
representations for time and energy. The obstruction is the need of com- 
pleteness of the g v in R and of the f v in R + to build systems of generators 
for square integrable functions. This is a necessary condition to arrive at a 
POV measure. If in addition we find orthogonality, then we would reach a 
PV measure. Notice however, that these measures would correspond to fuzzy 
distributions around a central value, in much the same way that the Husimi 
Gaussian wave packets do. This can be seen with our g v and f u , which are 
not eigenstates of time or energy. 

tg v (t) = i(3g v {t)- ig v -i(t), u f v (u) = v f v+1 (u) (17) 

The construction of eigenfunctions is hopeless at this stage. For instance, be 
ipioj) = J2v>v c u f u (u) a putative eigenfunction, so that = \ip{uj). 

Then ip{u)) oc 5{uj — A); expanding this in the f v {uj) would at least require 
that they form a complete set, which is not the case. The power of this 
simple example is that it signals the way to proceed. 



Energy and time representations 

Orthogonal polynomials are a useful tool to solve a variety of problems in 
physics. The Laguerre polynomials L°(x) constitute a set of orthogonal 
polynomials on the interval (0, oo) with weight x a e~ x . Accordingly, we can 
define the set of orthogonal functions 

= {0 if uj < 0, else c a n {—) a ' 2 e^ 2 " L a n (u/u )} (18) 

UJq 




= y j n(-l) m ( n + a \ (19) 
to \n-m) ml 

and ujq is a short of width of the exponential distribution that also serves 
to keep dimensions right (still, H — 1). This definition, provides a discrete 
denumerable set of modes that constitute a basis for all the functions be- 
longing to L 2 (0,oo). These functions comprise all functions with positive 
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frequencies (i.e. uo > 0) that can also be associated to probabilities (in fact, 
/ °°(ia;|<^"(ct;)| 2 = 1). The orthogonality and completeness relations read 
explicitly as 

roc 

/ dutfttu) tfttu) = 5 nn , and £ <H <p%(u>) = 5(u> - a/) (20) 
The physical information contained in the modes is readily obtained. In 
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Figure 1: The lowest modes for a = 2. uo is given in units of c^o and t — r 
in units of uOq 1 . The number of maxima for each mode is n + 1. The mode 
widths grow with increasing n. 

this representation, the Hamiltonian, in spite of its simple form H = uo is a 
non-diagonal matrix with elements given through 

n+1 

#|n> = E dln"o\™) (21) 

m=n— 1 
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In the same way 

n+2 

H 2 \n) = ]T d 2 mn ul\m) (22) 

m=n— 2 

The coefficients d m nS are computable by standard methods. Notice that d 2 mn 
is not (d^) 2 . The average value of the energy in the n-th mode is 

(n\H\n) = d\ n = (a + 2n + 1) oo (23) 

and the average of the squared Hamiltonian 

(n\H 2 \n) = d 2 nn = [(n + l)(a + n+ 1) + (a + 2n+ l) 2 + n(a + n)} u 2 (24) 

Finally, the variance of H, (AH)^ in the n-th mode is given by 

(AH) 2 n = (n\H 2 \n) - (n\H\n) 2 = d 2 nn - {d\ n f 

= [(n + l)(a + n + 1) +n(a + n)]u% (25) 

Another question to investigate are the eigenvalues E and eigenfunctions 
of the Hamiltonian. ^ E can be expanded in terms of the modes as ^>e — 
J2n c n | n ) which translates into the set of linear equations J2n d x mn c n = 
E c m left to the interested reader as an exercise. 

We take the next step forward and define wave packets centered at the 
point (r, E) in time energy space. To simplify the discussion we assume that 
E is the average energy in the n-th mode E — (a + 2n + 1)cj - We assume - 
for our purposes here - that we can trade (r, n) by (r, E). Finally, the wave 
packet we are looking for is: 

= e-») (26) 

The meaning of this time r just introduced can best be explored in the 
time representation. The Fourier transform of the modes \n, r) bring them 
from the energy to the time representation that is 

I fOO 

= -7j= / ^^" T V:H (27) 
V27T Jo 

Notice that ip% T (t) = i J no(^~ T ) by construction. Also, as the (p° T (iv) are 
a basis in R + they constitute a basis for L 2 (R). So, we conclude that they 
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form a system of imprimitivity leading to a PV measure. This is valid for 
both, the uo and the t representations as they are unitarily equivalent. 

The explicit expression of the modes in the t representation is quite cum- 
bersome, but throws some light on their physical content and on the way 
that the two different spectra are connected. From ref jU] we learn that the 
integral in (|27|) can be carried out giving: 

- "<> a T(a/2 + l)T(a + n + l) 
^ nA ) ~ n n!r(a + l) 

2 F l (-7i, a/2 + 1; a + 1; (1/2 + iu {t - r))' 1 ) 

{1/2 + iu Q {t-T)yi^ 1 ' 

Being the first entry a negative integer — n, the hypergeometric function 
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terminates. Explicitly 



2 Fi (-n, a/2 + 1; a + 1; (1/2 + ^ (t - r))" 1 ) = 

t ( 7 )m i n 2 + ^ (1/2 + ^ o(t " r)rm (29) 



where each term in the sum in ()29j) comes from the corresponding term in 
the sum in 

Due to the fact that ^"oW = "^nol - 0> we have (n, r|T|n, r) = r. To 
compute (n, r\T 2 \n, r), we consider the derivative of the mode (p"(uj)(uj = 1): 

£v» - (|j - i) v>» - I^^^Me^, o), (so) 

where we define 0(n, 0) = 1 when n > and zero otherwise. On the other 
hand, the Laguerre polynomials verify 

n 

l:=y: L a m -\ (si) 

m=0 

and thus 

d n 1 r a 

-r-<P» = E JWV) - r<(u;) - ^0)6(1, 0), (32) 
aio l=Q a c n _ 1 



where 



a c a 



Ki = l,^n-l + l). (33) 



2ci 



The matrix element (n,T\T 2 \n,T) results, after some straightforward opera- 
tions, 



oo roo 

2„/,q U\J.a fj.\ _ ,2 , / j++2„ 



(n,T\T 2 \n,r) = / dtt 2 ^ r (t)^ T (t) = r 2 + dtt 2 ^ (t)^ (t), (34) 

J — OO J —CO 

where we use (n, 0|T|n, 0) = 0. 

The last integral in (|34j) is given in energy representation by 



roo 



-oo 



d 2 _ r°° , r d 



dtt^ (t)^(t) = - 1 <H = y o du\— <(^)] 2 .(35) 
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Figure 3: (AT)„ (in units of cUq x ), as a function of a, for n = (continuous) 
and n = 3 (dotted). 

Making use of (jBSj) we get 

1 n / r a \ 2 

(n, r\T 2 \n, r) = r 2 + - + £(iV«) 2 + -M (n, 0) 

4 Z=0 \ C n-l/ 

- 2^_e(n,0)5:jvs/ d^rVK-iM (36) 

c n -i /=0 J0 

These expressions can be readily computed. The results confirm the expec- 
tations: The value of T 2 is near r 2 , the closer to it the larger a. We plot this 
behaviour in Fig. 3. Notice that (|3*K|) is nothing else but the variance (AT) 2 
for the mode n, something explicit after (|34jl . 

It is possible to build minimal uncertainty packets by combining differ- 
ent modes with appropriate coefficients. Instead, we will show how close 
the packets are to the lowest uncertainty. We do this in Fig. 4. Perhaps, 
the most prominent feature of these uncertainty relations is that the lowest 
uncertainty reached for each mode is bounded by n + 1/2, an asymptotic 
value corresponding to a — >• oo. On the other hand, the uncertainty remains 
bounded for the cases of physical interest. 

Finally, we look for the time of arrival at x of some state |\? ). We describe 
the state by a ket to avoid unwieldy notation, but the discussion could apply 
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Figure 4: (AH) n (AT) n as a function of a, for n = 0, 3. 



equally to density matrices. Working in parallel to what done to get and 
(JHJ), we could give the modes for mean time of arrival r at x with direction 
s = ±1 (for right and left movers respectively). They are given by: 



poo 

(nrxsl^) = / duj (x \p(oj) ) tp%Juj) (u I*) 

JO 



(37) 



the notation p(u;) stands for the dispersion relation governing the system. 
For a massive non relativistic particle uj = p 2 /2m, etc. We choose the case 
of a massless particle where p(u>) = soj. This not only avoids square roots 
but focus on the very interesting case of photons. In fact for a photon [TQ] 
we could use = F(x) = E(x) + iB(i), whose Schrodinger equation 

reads 

.dF{x,t) 



dt 



(a A V) F(x,t) 



(38) 



where a is the photon spin. Admitting the quantum leap involved in reducing 
this to one space dimension, it would translate into i^- = — i ^ . This is 
what we are using for the dispersion relation. By decomposing the state in 



terms of the modes (a; | \l/ ) = J2 



m Ymrm 



one gets after some computation 



12 



.08 



.06 




-30 -20 -10 10 20 30 

T — SX 



Figure 5: \(mTxs\n)\ 2 for m— 1, n = 3 and a = 2. 



that 



(n r x s | ^ ) 




{iu (r - sx)) 



m+n 



(1 + IUJq{t - sx)) 



m+n+1 



(39) 



-n; 



-m 



n 



a; 



1 + (a;o(r - sx))' 
(a> (r - sx)) 2 



where = c" 0^ r( -" 1 ~^, Q!+1 ' > . In Figs. 5 and 6 we plot these quanti- 
ties for two cases of interest. Notice the distribution of of these values around 
t = sx. This was expected. Notice also the zero at this value. This is just 
the modes orthogonality. 



Conclusions 



We have constructed a representation for time and energy in which both 
operate in a selfadjoint manner. The representation has the nice property of 
being suitable for probabilistic interpretation. It is given in terms of a set 
of modes build in terms of Laguerre polynomials and their weight functions. 
This allows to surmount quite easily the difficulties associated to the different 
support of the energy and time spectra. These difficulties not only translate 
into the Pauli theorem, but also prevent some asymptotic behaviour of the 
modes, precisely the more interesting from the physical point of view [TT] . 
According to the Paley Wiener theorem XII jT2j, it is not possible to get 
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2 for m — 1, n = 3 and a = 20. 



exponential asymptotic behaviour for the functions of time that are Fourier 
transforms of bounded functions of energy. Here, we dealt with the problem 
by using our modes. They are not exponential in time; this is out of reach. 
However, they have a controlled variance and uncertainty relations, in a form 
that makes them suitable for physics. 



This work was partially supported by the Ministerio de Education y Cien- 
cia of Spain under project BMF 2002-00834. The work of L. Lamata was 
supported by the FPU grant AP2003-0014. 

Note added in proof: After the submission of this paper the authors be- 
came aware of the interesting article Jose M. Isidro, Phys. Lett. A 334, 370 
(2005), in which a similar problem is studied from a different approach. 
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